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Abstract 

Parton equilibration for a thermally equilibrated but chemically 
non-equilibrated quark-gluon plasma is presented. Parton equili- 
bration is studied enforcing baryon number conservation. Process 
like quark - flavour interchanging is also considered. The degree 
of equilibration is studied comparatively for the various reactions / 
constraints that are being considered. 
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It is known [ 1 ] that the initial system produced at RHIC energies have a finite 
nonzero baryon number density. The anti particle to particle ratio at mid ra- 
pidity for y/SNN = ISOGeV Au-Au collisions at the STAR collaboration, BNL 
report a noticiable excess of baryons as compared to anti-baryons as reflected by 
the yields of p /p { hovering between 0.6 and 0.8 ), A /A ( hovering between 

0. 7 and 0.8 ) [ 2 ] for varying transverse momentum, centrality and rapidity. 
In a subsequent reporting of the BRAHMS Collaboration, BNL [ 3 ], it is seen 
that during Pb - Pb CoUisions at ^/snn = 200GeV , the p /p ratio is a 
maximum of about 0.8 at zero rapidity and falls off for higher rapidity values. 

This excess of baryons ( and hence quarks ) over anti-baryons ( and hence 
anti-quarks ) clearly necessiates inclusion of a chemical potential into the the- 
oretical framework . This shall be the main point of emphasis of the present 
piece of work. 

Inclusion of a chemical potential into the framework of parton equilibration 
studies has been rather recent. Although in 1986, Matsui et. al. [ 4 ] reported 
on strangeness equilibration rates at nonzero chemical potential, it was not until 
the very end of the past century [ 5,6 ] that it was explicitly used in studying 
chemical equilibration processes. A chemical potential for massless quarks that 
equalled the system temperature at all temperatures was assumed by these 
authors. This however was not very realistic as it did not take into account 
baryon number conservation. Furthermore, they did not distinguish between 
quarks and antiquarks. 

In 2004, He et. al. [ 7 ] undertook a much complete study including baryon 
number conservation. They used an expansion of the number densities in powers 
of the chemical potential. They started with a quark- antiquark distinction but 
that was soon put away. 

The present work is in the same line, but with the following points of differ- 
ence: 

1. We explicitly include the finite strange quark mass. 

2. We maintain the quark-antiquark distinction. 

3. We include the full phase space calculations for the pair production and 
flavour changing processes in contrast with the factorised rates of many earlier 
works. [ 5,6,7,8,9,10 ] 

4. For a given initial baryon number density we iterate to obtain the initial 
value of the quark chemical potential which again evolves obeying the baryon 
number conservation equation. 

5. We study comparitively the degree of equilibration achieved for varying ini- 
tial conditions. 

At this stage, let us note the following points as regard to the chemical po- 
tential: 

1. Strangeness being very nearly conserved in strong interactions and the ini- 
tial strangeness content of the QGP flreball being zero it implies that strange 
quark and antiquark are always produced in pairs. This clearly indicates that 
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the strange quark chemical potential is zero. 

2. Considering light quark-antiquark pair production processes we can argue 
that since the net chemical potential on either side of a chemical reaction ought 
to be the same and gluon chemical potential is zero, this would necessarily imply 

f^q = -fJ'q 



Again, since the mass of u and d quarks are nearly equal and much less than 
that of the strange quark, we can safely treat the Hghter flavours in the same 
footing without any appreciable error being introduced. By a suffix q we shall 
indicate the lighter quarks in general and treat them as massless. 



1 Basic Thermodynamics of the System 

The distribution functions for the constituent partons of the chemically non- 
equilibrated state are taken to be as 

■^3 = ^elT_ Y ( ^ 

_ ^q{q)^ ^ ( lb ) and 

•'9(9) ge/T j^l ^ ' 

■fs^h^ ^,,t\^ ( 10 

where the notations have usual meanings. 

Of these, the first is the usual Bose distribution function [ 8 ], the second is 
the Modified Fermi-Dirac type distribution function with an exponential term 
in the numerator to include the light quark chemical potential [5 ] and the third 
is the usual modified Fermi-Dirac distribution function for the massive strange 
quark [ 8 ] , each being scaled by the non-equilibrium fugacity in order to describe 
the degree of equilibrium achieved for a chemically equilibrating system. 

Using standard rules of Statistical Mechanics [ 11 ] we can obtain the number 
density n , energy density e and pressure p of the system. The results are 
generally represented as : 

t — tg + n f(tq + tq) + 2ts ( 2 ) where t = n,e,p with 

.16C(3), 9C(3)n/,^ ii, ^ 



^^^D-D'-W^' (-) 



1-i.q 



£ = i^^a + ^^{\e- + Xqc -)+ 

AA..^f:(-l)-H^f^ + ^^}]T^ (2b) 

71-2 {kXsY (kxs) 
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\ I / \ ^ T 



(.c) 

where gives the number of massless quark flavours (=2 for our case) and 

The baryon number density equals ( for two massless quark flavours) 

rif, = 2X — {nq — n-q) ( 3 )• For given initial values of baryon number 
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density, non-equilibrium fugacities and temperature we can iterate to obtain 
the initial value of the chemical potential. The baryon number conservation 
equation is given by d^j,{nbu^) = (4) 

which may be expanded to obtain the rate of change of the chemical potential 
with time. We obtain 

^iq = \Bi + XqB2 + TB^ + Bi (5 )where 

Bi = 77= ( 5a ) 

B2 = n /T , \ -u IT ( ) 

Xqef^"'^ + XqC 1^1' ^ 

B3 = ^ - i " ,^ " ^ ( 5c ) and 

nXqe^J^^X-qe-^^/^) 

" ^(v^v^Tv^^ ^ ^ 

It is to be noted here that for a baryon-free plasma ^iq vanishes. As we shall 
see shortly, due to these non-zero coefficients, the Hght quark and anti-quark 
number density evolution equations get coupled with one another. 



2 Parton Equilibration Equations 

The fundamental equation that dictates the parton number density evolution is 
given by 

0^[nkU^') = + — = [Rgain - Rloss) 

(6 ) 

where the RHS gives the difference of rate of gain and loss of the parton species 
k i.e the net rate of change of the number density. Using standard procedures 
[ 10, 11 ] this equation leads to the parton number density evolution equations 
as we shall see shortly. 



4 



2.1 Gluon Number Density Evolution Equation 



The Gluon Number Density Evolution Equation is given as 

d^ingU^') = [Rgg^ggg - Rggg^gg) - Y^iRgg^U " ^ii-^gg) 

i 

(7 ) 

where the sum is over all quark flavours present. Substituting for the gluon 
number density and substituting for the rates we have 

AgGi+TGs + G3 = 

(8 ) 

where 

Gi = l/Xg (8a ) 

G2 = 2,/T (8b ) and 

G3 - (1/r) ~ (i?3(l - \g) ~ E^ R2^/ng ) (8c ) 

where, as in earlier works [ 8, 10, 11 ] we have introduced rates i?2i and R3 . 
We shall look more closely at these rates soon. 



2.2 Quark and Anti-Quark Number Density Evolution Equations 



As mentioned earlier, due to the coefficients Bi , the light quark and anti- 
quark number density evolution equations get coupled to each other. We obtain 
the following results: 



2.2.1 Massless Quarks 



The equation is given by 

XqQi+XqQ2 + tQ3 + g4 = (9 ) 

where 

Qi = 1^ + ^ ( 9a ) 

Q2^^ (9b ) 

Q3 = |; +^-^ (9c) 
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Q4=- + ^- SQ (9d )with 

T 1 

SQ = {{Rgg^qq- Rqq^gg)/nq] - {{Rqq^s-s - Rss^qq)/nq} (9e ) 

for light quarks 



2.2.2 Massless Anti- Quarks 



For the anti-quark case, we obtain the equation as 

\gAQi+\gAQ2 + TAQy. + AQ^ = ( 10 )where 

AQi = -^ (10a) 

1 B7 

= 3^ - ^ ( ) 

AQ, = ^ + ^ (10c) 

AQ4 =--^- SaQ ( lOd )with 

T T 

SaQ = {{Rgg^qq - Rqq^gg)/nq} - {{Rqq^ss " RsS-^qq) / Tlq} ( lOe ) 

for light anti-quarks 



2.2.3 Massive Strange Quarks 



For massive strange quark the equation is given by 

XsSl+tS2+S3=0 ( 11 ) 

where 
Si = ^ ( 11a ) 

3t[ ikx^y (k^s) ^ 



( lib 



S3 =- - SQs ( 11c ) with 

T 

SQs^ {{Rgg^ss - Rss^gg)/ns} + 2{{Rqq_s-j- Rss-,qq)/ns} (lid) 

for strange quarks 
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2.3 Energy-Momentum Conservation Equation 



From the energy- momentum conservation equation d^T^'^ ~ ( 12 ) 

we obtain for (1+1)D — + = ( 13 ) 

dr T 

( which we shall refer to as Bjorken's equation henceforth). We get on substi- 
tuting for energy and momentum 

Th + \jh + \h + V9 + A./s + /lo = ( 14 ) 

where 



72 4 ,K2{kxs) 5Ki{kxs) Ko{kx 



-m 



^2T""^'^^^ {kx^Y + 12(fcx,) + 12 ^ ^ ^^^^ 

/7 = ( 14b ) 

/. = ;^™lELi(~l)'=-H^t^ + ( 14c) 



/s= ^e^T^ +^(A,e^-A,e-^)4iT^ ( 14d ) 

/9= ^e-^T^ +^(A,e^-A,e-^)tr4 ( 14e 



Ao = l^A,T^ + ^(V^+A,e-^)T^- 
Swf, v°<= ^ ^^k-lf ^K2{kxi) ^ Kijkxi) ^ 

+ ?;^(A,e^-A,e-^)^T^ (14f ) 



There remains a very important point to mention regarding the evaluation of the 
rate of change of the chemical potential. Since the baryon number conservation 
equation is obtained directly from the quark and antiquark number densities, 
as are the massless quark and anti-quark number density evolution equations, 
these two evolution equations (coupled by the four B - coefficients) alongwith- 
with the energy-momentum conservation equation ( with substitutions for As , 
Ag from the respective number density evolution equations) do not gener- 
ate a set of solvable coupled equations for A, , A^ , T as the system 
deeterminant vanishes. 

To tackle this problem, we adopt an approximation scheme. While evaluat- 
ing Bj, [j = 1, 2) we use a truncated form of the expansion of the exponential 
in the numerator keeping the rest of the expressions unchanged. 
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3 Parton Equilibration Rates 

3.1 Gluon Multiplication Rate -R3 

The gluon multiplication rate has been calculated by Xiong et. al. [ 12 ]. By 
explicitly calculating the matrix element [ 13 ] (summed over all the final states 
and averaged over all initial states) we can obtain the gluon multiplication rate. 
However, to avoid the huge calculations of evaluating 25 Feynman diagrams [ 
13 ] involved, we fall back on the treatments used in earlier works [ 6,11 ]. We 
postulate that the gluon multiplication rate depends on the chemical potential 
via the Debye Screening mass. 

The Debye Screening mass suitable for a multicomponent chemically non-equihbrated 
parton plasma is given by [ 14 ] 



'D^%J dkk[NJg + j:j,] ( 15 



where the sum runs over all flavours i, while Nc gives the number of colours. 
To accomodate for antiquarks and remembering that our number of flavours is 
3 and not 6, we propose the following modiflcation: 



m 



D 



■J dkk[3f,+ 1 E.=„..,.(/. + /-)] (15a) 



Using standard techniques [ 5,6,8,11 ] we can obtain the following result for the 
mean free path A/ : 



16) 

which for zero chemical potential using = AnasT'^Xg ( 17 ) 

reduces to the well known [ 8,11 ] result: 

'f' = ^^"-^ 1 + 8.1a,/9 ( 1« ) 

Using standard methods [ 5,6,8,11 ] we get the modifled differential cross section 
as 

- w - < " ' 

RecaUing the deflnition of i?3 = -^(^sng ( 20 ) with 

(J3 ~< (7{gg ggg)v > ( 21 ), the thermally averaged velocity weighted 

cross section, we obtain the required rate as 

AK/(A,) ( 22 ) 



i?3 27a3 



T 2 
where 

s 

I(Xg) = / dx dz- -\ ^ = + 

Ji Jo 1 + zy x^[x + (1 + z)xd]^ - 4.XZXD 



cosh ^y/x 



with 



^[l + x(l + 2)yD]2-4xzyi5 



(23) 
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(24) 



3.2 Quark Anti-Quark pair production rate 



We have parton production rates in the RHS of the number density evolution 
equations. Let us evaluate the quark-antiquark pair production reaction rate 
i?2g . We have [ 4 ] 

R2,^RlL-RiL (25) 

where 

' J {2^f2Ej i2nf2E2 J {2ttY2E:, J {27rf2Ei 

) 

and 



( 26 



' {2nf2Ej {2TTf2E2 J {2TTf2E3 J (2^)^2S4 
S\pi+P2-P3-P4)nMgg-./i^ + faiPiW + h{P2))IMk{Pi) - ( 25b 

) 

Following [ 4 ], we can say that there are three topologically distinct Feynmann 
diagrams that contribute towards the quark-antiquark pair production process. 
Evaluating them performing traces and finally adding them up we can find the 
net squared matrix element. We basically follow the Hnes of [ 4 ]. Transforming 
variables as 
q = Pi+P2 

P = i^iPi -P2) 
9' = P3 + Pa 

P' = ^iP3-Pi) 

with restrictions 

qo > 2mi 

s = 9o - |q| > 4mf 

Art 4(1 -i^) 

and transforming the three dimensional integrals to four dimensional integrals 
using 

d^Pi 
2E, 

with the new set of variables 



J jdS.6(p'^m^) (28) 
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qo = —Tlnv + 2mi 

9.-1 4m? 1 ( 29) 

2, q' 

we arrive at the rate ^ ^ ^ 2 2 

Rlg^-^Tj du j dv I dx I dy (1 -)Hqo - 4.171^)^ f Quarks fphasel 

Jo JQ Jo Jo ^' S 

(30) 
where 

fQuarks = /g(y +po)/s(y -Po)(l - /q(y +Po))(l - My -Po))" 

(1 + /,(f +Po))(l + /,(f -Po))/,(f +p(,)Mf -Pf,) ( 30a ) 

and ^ A A 

fphasel =A + + + ^ij^ + ( 30b ) 

with 

A = 3[l-[l-M][(i^^!fc^+.V]]-|-24^ (30c) 

B = f [1+^-^ + 4] ■-■-(30d) 
%0e) 

i^^ = [l_[l_M][(i„^2_^2)^2[l-^]ia;j/]i ( 30f ) 

s s 
4m? ,1 , , 

A± = 1± [1 '-]2xy ( 3Qg) 



3.3 Quark Flavour Changing Rate 



For the quark flavour changing process we have [ 4 ] 
where 

5^ij)l + P2 - P3 ~ Pi)AM,-^^qq\ /g(Pl)/g(P2)(l " /s(P3))(l " fs{Pi))- ( 31a 

) 

and 

Rf = / / / / -^i2.f 

J {2nf2EiJ {2TTy2E2 J {2TTf2E3 J (2^)'^2S4 

SHpi+P2-p3~P4)nMsl^gg\\l ^ fgiPlW - fq{p2))fs{p3)f-s{Pi) - ( 31b 

) 

Following [ 4 ], we can say that there is only one type of topologically distinct 
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Feynmann diagram that contributes towards the quark flavour changing process 
or the strange quark pair production process. Evaluating it, performing trace 
calculations we can find the squared matrix element. We basically follow the 
lines of [ 4 ] . Transforming variables as in the case before and performing iden- 
tical operations we can get the rate as 

^qg^TT^^ / / dy—{l {Qq - 4:mf)'^ f strange fpha8e2 

'^'^ Jo Jo Jo Jo ^ s 

(32) 

with 9o 9o 9 1 

fQuarks = .^(y + Po)/q (y " Po) (1 - ./s(y +Po))(l ~ My-p'o))- 

(l + /.(f +Po))(l + /.(f -Po))/,(f +P^,)Mf -P^,) ... (32a 

) 

u.as.. = - + .^,^]] + M ...( 32b ) 

s 2 s 



4 Results 

4.1 Initial Conditions 



As an initial condition at the point of thermalization ( i.e the point after which 
the system evolves according to the laws of hydrodynamics) we have predictions 
from the HIJING and SSPC models [ 15 ]. In addition to these inputs we 
shall require inputs for at least two of the following: initial baryon number 
density, initial light quark chemical potential and ratio of initial non-equilibrium 
fugacities of light quarks and antiquarks. 

In absence of initial values thereof ( although we have some initial conditions 
calculated at n = 0.1 fm/c [1,16 ], which is at a much earher time than the 
time of thermalization at Ti = 0.25 fm/c for SSPC initial conditions and 
earher still for HIJING at r, = 0.6 - 0.7 fm/c [ 15 ] ) let us study the 
relative degree of equilibration for varying initial conditions although it must be 
emphasized that most of these initial conditions would be of purely academic 
interest only, as they cannot be realised in practice at the collider experiments. 
Nevertheless, it is hoped that this exercise would pave the way for a clearer 
understanding of the physical processes involved. 

Even though two inputs remain rather arbitrary, we use both set of initial 
conditions ( LHC and RHIC ) to study the trends. For the usual inputs we use 
the following sets of data ( SSPC ) 



LHC 



n{fm./c) T{GeV) Xg 




A. 


0.25 1.02 0.43 


0.086 


0.043 



RHIC 
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Tiifm/c) 


T{GeV) 


^9 


A, 


A. 


0.25 0.668 0.34 


0.064 


0.032 



As a representative plot we use baryon density 0.15/fm"3 and ratio=1.5 with 
the initial conditions as given above. We also try to see the effect of including 
the quark flavour changing processes ( qfcp) . 

4.2 Observed trends 

We study comparatively the outputs with given initial baryon number density 
and initial light quark to anti-quark fugacity ratio. We vary the baryon number 
density between 0.21/fm"3 and 0.11/fm"3 while we change the ratio between 
1.9 and 1.1. We arrive at the following conclusions : 

1. The nature of variations of the physical quantities are more or less in 
the line of earlier works. The results for RHIC and LHC initial conditions are 
shown in figures 1 and 2. Temperature, Non equihbrium fugacity and Chemical 
Potential variations with RHIC initial conditions are given in figure 1 while 
that for LHC are given in figure 2. The decaying curves give temperature 
variations while the positive rising curves are sequentially (from top) for gluon, 
light quark, light anti-quark and strange quark fugacity variations. Plots for 
both quark flavour changing process included and excluded cases are shown. 
For identification please see point 4 below. 

We observe that 

i) As expected, the temperature falls with time while the non-equilibrium fu- 
gacities increase. 

ii) Contrary to the Juttner case, the chemical potential remains negative all 
along and as expected approaches zero as the system equilibrates. 

iii) The QGP,as expected remains to be gluon dominated. 

iv) For a system at higher chemical potential, the temperature falls at a slower 
rate signifying lesser amount of energy expenditure to create partons which 
shows up in the slower rise of all partons except for the light anti-quark, which 
shows a larger growth rate due to the presence of the exponentiated chemical 
potential. 

2. For a given initial ratio, except for the light quarks and anti-quarks, the 
output does not depend much on the initial baryon number density. For the 
Hght quarks and antiquarks this variation is due to the presence of the respective 
non-zero chemical potentials. 

3. For a given initial baryon number density we can recast the equation for 
baryon number density in the form 

„ c n 

pfiq/-i _ 1 / 1 

2A, V ^ 4A2 



12 



( 33 ) 

where C is a constant for a given temperature, light quark fugacity and given 
baryon number density. Here D is the light quark to antiquark initial fugacity 
ratio. For a fixed Hght quark fugacity and temperature, as D falls clearly the 
RHS of the above equation increases which indicates a rise in the chemical po- 
tential. Again for a system of higher chemical potential, the temperature has 
to drop at a slower rate due to the constraint imposed. Hence we observe that 

i) The temperature falls at a slower rate for a smaller value of the light quark 
to antiquark initial fugacity ratio. 

ii) As the temperature falls at a slower rate, it would imply a lesser expenditure 
of energy to produce particles in general . This would show up in the slower 
rise of all fugacity values except the light antiquark, which would show a higher 
growth rate. This is due to the exponentiated chemical potential part. 

4. For inclusion of the quark-flavour-changing process we observe the fol- 
lowing: 

Due to the additional production of s-quarks from massless quarks via qfcp, 
we see an additional increase in fugacity of the strange quark while, the rate of 
equilibration falls for all other non-strange quarks . 

5. For the chemical potential of the light quarks, clearly the value thereof 
tends to level off as the system equilibrates . The chemical potential remains 
negative all along, which can be attributed to the choice of the distribution 
function. 
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Figurel: RHIC Case 
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